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We discuss the phenomena of symmetry non-restoration and inverse symmetry breaking in the 
context of multi-scalar field theories at finite temperatures and present its consequences for the 
relativistic Higgs-Kibble multi-field sector as well as for a nonrelativistic model of hard core spheres. 
For relativistic scalar field models, it has been shown previously that temperature effects on the 
couplings do not alter, qualitatively, the phase transition pattern. Here, we show that for the 
nonrelativistic analogue of these models inverse symmetry breaking, as well as symmetry non- 
restoration, cannot take place, at high temperatures, when the temperature dependence of the 
two-body couplings is considered. However, the temperature behavior in the nonrelativistic models 
■ allows for the appearance of reentrant phases. 
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The study of symmetry breaking (SB) and symmetry restoration (SR) mechanisms have proved to be extremely 
useful in the analysis of phenomena related to phase transitions in almost all branches of physics. Some topics of 
current interest which make extensive use of SB/SR mechanisms are topological defects formation in cosmology, the 
£f) ' Higgs-Kibble mechanism in the standard model of elementary particles and the Bose-Einstein condensation (BEC) in 
condensed matter physics. An almost general rule that arises from those studies is that a symmetry which is broken at 
zero temperature should get restored as the temperature increases. Examples range from the traditional ferromagnet 
to the more up to date chiral symmetry breaking/restoration in QCD, with the transition pattern being the simplest 
one of going from the broken phase to the symmetric one as temperature goes from below to above some critical value 
and vice-versa. 

However, a counter-intuitive example may happen in multi-field models, as first noticed by Weinberg |l( who 
considered an O(N^) x 0(N^) invariant relativistic model with two types of scalar fields (with N$ and N$ components) 
Oh, and different types of self and crossed interactions. Using the one-loop approximation he has shown that it is possible 
for the crossed coupling constant to be negative, while the model is still bounded from below, leading, for some 
parameter values, to an enhanced symmetry breaking effect at high temperatures. This would predict that a symmetry 
which is broken at T = may not get restored at high temperatures, a phenomenon known as symmetry non restoration 
(SNR), or, in the opposite case, a symmetry that is unbroken at T — would become broken at high temperatures, 
thus characterizing inverse symmetry breaking (ISB). Here, one could argue that SNR/ISB are perhaps just artifacts 
of the simple one-loop perturbative approximation and that the consideration of higher order terms and effects like 
the temperature dependence of the couplings could change the situation. To answer this question the model has been 
re-investigated by many other authors using a variety of different methods with most results giving further support to 
the idea (see, e.g., Ref. Q for a short review of SNR/ISB). For example, the SNR/ISB phenomena were studied using 
the Wilson Renormalization Group |3j and the explicit running of the (temperature dependent) coupling constants 
has been taken into account, showing that in fact the strength of all couplings increase in approximately the same 
way with the temperature. This analysis shows that once the couplings are set, at some (temperature) scale, such as 
to make SNR/ISB possible, the situation cannot be reversed at higher temperatures. Two of the present authors have 
also treated the problem nonperturbatively taking full account of the cumbersome two-loop contributions The 
results obtained in Ref. Q were shown to be in good agreement with those obtained using the renormalization group 
approach of Ref. [3j and, therefore, also support the possibility of SNR/ISB occurring in relativistic multi-scalar field 
models even at extremely high temperatures, where standard perturbation theory would break down. 
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The mechanisms of SNR/ISB have found a variety of applications. For instance, in cosmology, where they have 
been implemented in realistic models, their consequences have been explored in connection with high temperature 
phase transitions in the early Universe, with applications covering problems involving CP violation and baryogenesis, 
topological defect formation, inflation, etc For example, the Kibble- Higgs sector of a SU(5) grand unified theory 

can be mimicked by considering the case N$ — 90 and N$ — 24 and has been used to treat the monopole problem 
0,EJ 13 • Setting — = 1 the model becomes invariant under the discrete transformation Z2 x Z2- The latter 
version has been used in connection with the domain wall problem J| . Most applications are listed in Ref. which 
gives an introduction to the subject discussing other contexts in which SNR/ISB can take place in connection with 
cosmology and condensed matter physics. These interesting results from finite temperature quantum field theory raise 
important questions regarding their possible manifestation in condensed matter systems which can be described by 
means of nonrelativistic scalar field theories in the framework of the phenomenology of Ginzburg-Landau potentials, 
like, for example, in homogeneous Bose gases Q. As far as these systems are concerned, we are unaware of any 
applications or studies of analogue SNR/ISB phenomena in the context of nonrelativistic scalar field models. 

In the context of condensed matter physics more exotic transitions are well known to be possible and similar 
phenomena to ISB/SNR have been observed in a large variety of materials. One of the best known examples is 
the symmetry pattern observed in potassium sodium tartrate tetrahydrate, KNa(C<tH^O§)AH20, most commonly 
known as the Rochelle salt, which goes |icj, as the temperature increases, from a more symmetric orthorhombic 
crystalline structure to a less symmetric monoclinic structure at T ~ 255A". It then returns to be orthorhombic phase 
at T ~ 297A, till it melts at T ~ 348A. It thus exhibits an intermediary inverse symmetry breaking like phenomenon 
through a reentrant phase. Other materials which arose great interest recently due to their potential applications 
include, for example, the liquid crystals [ll] and spin glass materials which exhibit analogous phenomena of 
having less symmetric phases at intermediary temperature ranges, known as nematic to smectic phases (ferro and 
antiferro electric and magnetic like phases), and compounds known as the manganites, e.g. (Pr, Ca, Sr)MnOs, which 
can exhibit ferromagnetic like reentrant phases above the Curie (critical) temperature |13 j . Actually, in the condensed 
matter literature we can find many other examples of physical materials exhibiting analogue phenomena of SNR/ISB. 
This same trend of the emergence of reentrant phases also seems to include low dimensional systems 01 • A discussion 
on these inverse like symmetry breaking phenomena in condensed matter systems has been recently summarized in 
Ref. |15| . Another motivation for the present work is the growing interest in investigating parallels between symm etry 
breaking in particle physics (Cosmology) and condensed matter physics (the Laboratory) as discussed by Rivers [l6| 
in a recent review related to the COSLAB programme. One of the most exciting aspects of such investigation is due 
to the fact that condensed matter allows for experiments which can, in principle, test models and/or methods used 
in Cosmology. 

Here, our aim is to analyze a nonrelativistic model composed of two different types of multi component fields. 
To investigate SNR/ISB we consider a model possessing an U(l) x U(l) global symmetry, that is analogous to the 
0(N ( / ) )xO(N^) relativistic model studied in f° r ^<t> = = 2, including both one and two-body interactions in 

the potential. Further, by disregarding the bosonic internal degrees of freedom, the model is considered as representing 
a system of hard core spheres. In the analysis that follows in the next sections we do not claim that this simplified 
model described in terms of scalar fields with local interactions will be simulating the phases behavior of any of 
the condensed matter system cited in the previous paragraph, but just that it suffices, as a toy model, to show the 
generality of the possibility of emergence of reentrant behavior in some simple condensed matter systems which can 
be modeled by coupled multi-scalar field models. The chosen non relativistic model is also simple enough to show the 
differences and analogies regarding the phenomena of SNR/ISB which occurs on its relativistic counterpart. 

We will show that, like in the relativistic case, SNR/ISB can take place when thermal effects on the couplings 
are neglected. We then consider these thermal effects by computing the first one-loop contributions to the couplings 
finding that, contrary to the relativistic case, SNR/ISB cannot persist indefinitely at higher temperatures when all 
symmetries are restored. In summary, the possible phase transition patterns seem to be completely different for the 
relativistic and nonrelativistic cases when the important thermal effects on the couplings are taken into account. This 
paper is divided as follows. In Sec. II we review the original relativistic prototype model. In Sec. Ill we present a 
similar nonrelativistic model of hard core spheres with quadratic and quartic interactions. We show how SNR/ISB 
cannot occur for such a system when the temperature effects on the couplings are considered, but they can only 
manifest through reentrant like phases, with symmetry restoration always happening at high enough temperatures. 
Our conclusions and final remarks are presented in Sec. IV. An appendix is included to show some technical details 
of the calculations. 
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II. THE EMERGENCE OF SNR/ISB PHENOMENA IN THE RELATIVISTIC MODEL 



At finite temperature the relativistic multi-scalar field theory was first studied by Weinberg [l| who found evidence 
of SNR/ISB taking place at finite temperatures. On his work, he considered a prototype model composed of two 
types of scalar fields, (f> and i/j with Na, and Nj, components, respectively, which is invariant under the 0(N<p) x O(iV^) 
transformation. Such a model has a lagrangian density which can then be written as 



= ~ -f<P - f (0 2 ) 2 + ^) 2 - -f <A 2 - ^0/> 2 ) 2 - • (2.1) 

The self-coupling constants A^ and and the cross coupling A in Eq. 12.1(1 are traditionally considered as all positive. 
However, it is still possible to consider A negative in 1(2.1(1 provided the potential is kept bounded from below. It is 
easily seen in this case that the boundness condition for the model 12.1(1 requires that the couplings satisfy 



A > 0, \^ > 0, A^Av, > 9A 2 . (2.2) 

The fact that the cross coupling, A, is allowed to be negative has interesting consequences as is seen from the 
one-loop thermal mass evaluation. As usual, the temperature effects on the zero temperature mass parameters mf 
(where i = <f> or -0) can be computed from the (thermal) self-energy corrections Sj(T) from which the thermal masses, 
Mf(T) = nij(0) + 'Ei(T) are obtained. The thermal masses have been first calculated with the one-loop approximation 
which, using the usual rules of finite temperature quantum field theory (see e.g. [H.ll7|') and in the high temperature 
approximation, m^/T, m^/T <C 1, leads to the results 



m2(t) 



12 



(2.3) 



and 



Ml{T) 



12 



A, 



1 /AL + 2 



X- 



(2.4) 



where we kept only the leading order relevant thermal contributions in the high temperature expansion of Ej(T), 
which will be enough for the analysis that follows. Note also that the zero temperature quantum corrections to both 
masses and coupling constants are divergent quantities and so require renormalization. This is done the standard 
way by adding the appropriate counterterms of renormalization in 1(2.1(1 (see also 3). We are only interested in the 
thermal quantities (that are finite) since the zero temperature quantum corrections to masses and couplings can be 
regarded as negligible as compared to the finite temperature contributions. In Eqs. ((2.3(1 . 1(2.4(1 as well as in the 
relations below, the mass parameters and and couplings A^, A</> and A are just to be interpreted here as the 
renormalized quantities instead of the bare ones. It is obvious, from the potential term in the lagrangian density ((2.1(1 , 
that if one of the mass parameters mf is negative the 0{Ni) symmetry related to that sector is broken at T = 0: 
O(Ni) — ► 0(Ni — 1). Thermal effects tend to restore that symmetry at a certain critical temperature, upon using 
Eqs. 1(2.3(1 and 1(2.4(1 . given by 



-12m- 



1 fN t + 2 



^ N 3 



1/2 



(2.5) 



However, if mf < 0, Eq. ((2.5(1 shows that for a negative cross-coupling constant, A < 0, and for |A| > Aj(A 7 , + 2)/ (3Nj), 
T c ^ cannot be real. In other words, the broken symmetry is never restored (SNR). At the same time if mj > 
(unbroken 0(N,) symmetry at T = 0), but A < and |A| > X i (N l + 2)/(3A,) then from Eqs. l(2T3|). l(2~4l and (|2~5|l . 
we can predict that, as the temperature is increased, the symmetry will be broken at T c , instead of being restored 
(ISB). For example, let us suppose that A < and 



(2.6) 



4 



In this case the boundness condition assures that |A| < \^{N^ + 2)/(3N$). Then, if < one has broken O(N^) 
symmetry at T = 0, but M|(T) will eventually become positive at the corresponding T c ^, given by Eq. Q2.5JI . 
restoring the symmetry. If > 0, then M|(T) > for all values of T and the model is always symmetric under 
0(Nip). On the other hand, if < 0, our choice of parameters predicts that the O(N^) symmetry is broken at 
T = and that it does not get restored at high temperatures, a clear manifestation of SNR. At the same time, if 
mk > 0, the 0(N t p) symmetry, which is unbroken at T — 0, becomes broken at a T > T c ^, which is a manifestation of 
ISB. Obviously, which field will suffer SNR or ISB depends on our initial arbitrary choice of parameter values. Note 
that when A = the theory decouples and SNR/ISB cannot take place. In this case one observes the usual SR which 
happens in the simple O(N) scalar model. 

An issue that arises, concerning the results discussed above, is that the coupling constants are scale dependent 
in accordance with the renormalization group equations. Therefore, at high temperatures not only the masses get 
dressed by thermal corrections but also the coupling constants, so we must answer whether the intriguing phase 
transitions patterns discussed above, for A < 0, can hold in terms of the equivalent running coupling constants. This 
issue was analyzed by Roos ; 3j, who used the Wilson Renormalization Group (WRG) to evaluate the Xi(T) and A(T). 
His calculations revealed that the strength of all couplings increase, at high T, in a way which excludes the possibility 
of SR in cases where SNR/ISB happen. He also showed that the running of coupling constants with temperature 
as predicted by the one-loop approximation, as adopted in the present work, is robust up to very large scales. In 
addition to that, the two-loop nonperturbative calculations performed in Ref. also support, from a qualitative 
point of view, Weinberg's one-loop results. 

As one notices from the equations which describe the thermal masses, Eqs. (H2J) and the appearance of 

SNR/ISB is directly related to the relation among the different couplings when A < 0. It is then useful to define the 
quantity 



which takes those relations into account. Then, in terms of temperature independent couplings, the critical tempera- 
ture, Eq. H2.5JI . can be written as 



-12M?^ 1/2 



One can easily see that SNR/ISB may occur when one 1 of the Aj is negative. 




FIG. 1: The diagrammatic representation for the effective coupling \#(T) up to the one-loop level. The continuous lines stand 
for the 4> propagators, while the dashed lines represent tp. 



One of the main results of Ref. Q states that SNR/ISB can occur in both sectors, for some parameter values, a situation which is not 
allowed at the one-loop level. 



FIG. 2: The diagrammatic representation for the effective coupling Ay,(T) up to the one- loop level. 



Let us now check the robustness of SNR/ISB when the effective, temperature dependent couplings are considered. 
The thermal effects on all the three couplings, at the one-loop order, are considered in terms of the corrections to the 
four-point 1PI Green's functions. All diagrams at the one-loop level contributing to the effective couplings X^(T), 

These diagrams, with zero external momenta 2 , are easily 



A^,(T) and A(T) are shown in Figs. Q]|21and|21 respective! 
computed at finite temperature (see for instance Refs. 
leading order, one obtains 



sly. ±J 

urn 



). Using again the high temperature approximation at 



8 A 2 , N i> X 2 
A + — A 



(2.9) 



A^,(T) ~ A,/, + 
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In 
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(2.10) 



and 



A 



A(T)^A+ — 



A + In ( — 

8vr 2 V M 



N,i, 



2\ x 1 

A.+ 2 



A i 



-M- 

Air 2 \M 



(2.11) 



where Mq is a regularization scale. In writing the above equations we are once again assuming that the tree-level 
couplings in Eqs. (|2.9I) . (|2.1U|) and l|2.11|l are the renormalized ones and we also are only showing the relevant 
high temperature corrections. The same expressions were also obtained by Roos in Q (note however that different 



normalizations for the tree- level potential as well as Mq = To, N<h 



Njp = 1 were used in that reference). In 0, 
the numerical solution of the one-loop Wilson renormalization group equations was also compared to the usual flow 
equations for the constants obtained from the one-loop beta-functions and shown to agree well with each other up to 
very high scales. The flow equations referring to the perturbative effective coupling constants Eqs. I|2.9|l . i|2.1U[l and 
(|2.11l) are expressed in term of the dimensionless scale T/Mq 01 as 



d\<f,(t) 
dt 



3 
8^ 



(2.12) 



dt 



3 
8^ 



(2.13) 



2 Recall that those are contributions to the effective potential, which generates all 1PI Green's function with zero external momenta. 
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where t — ln(T/Mo) was used. The solutions of the flow equations i|2.12[l . (|2.13|) and (|2.14|) . with initial conditions 
given by the renormalized tree-level coupling constants, can also be easily seen to be equivalent to the solutions for 
the set of linear coupled equations, 
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A(T) 



3 , ( T 
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N d> + ; 
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16tt 2 V M, 
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1 / A^ + 2 
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N* + 2 



A^A^T) - 
A^A^,(T) 



AA(T) 



AT, 



A (T) + 



2 



AA(T) 
2 



A,* + 



1 fN, 



A, 



4tt 2 VMq 



AA(T) . 



(2.15) 



The results obtained from the flow equations given above, or equivalently from the solutions of coupled set of equations 
(|2.15|) , are standard ways of nonperturbatively resumming the leading order corrections (in this case the leading log 
temperature dependent corrections) to the coupling constants. For instance, Eq. (|2.15() is exactly the analogous 
procedure used for the one-field case for summing all ladder (1-loop or bubble) contributions to the effective coupling 
constant. For the multi-field case, the perturbative approximation for (|2.15l) is again given by Eqs. (|2.9I) . (|2.1UI) and 
(|2.11l) at the one-loop level. In our case, these equations are useful to test how robust is the phenomena of SNR/ISB 
and will be used below in our analysis. Later, in the next section for the nonrelativistic limit of Eq. (|2.1|l . we will 
also construct the analogous of these nonperturbative equations for the temperature dependent effective couplings. 



In terms of the effective temperature dependent couplings. 
(|2.7|l becomes 



\<f>(T), X^(T) and A(T) the quantity analogous to Eq. 



A i (T)=X i (T)-l-^-\+\(T)^-, (2.16) 



or, more explicitly, using Eqs. I|2.9|l . (|2.10() and l|2.11[) . 



N z + 2\ , X JV,- , . 2 (N i + 8)(N i +2) 1 fT\ (N i + 2)X i + (N j +2)X j ^fT 



A i( T) = X, ( ^- ] + A f + \f - M „ I ) + AiV, ' — ' - i„ 



It is clear from the expressions for the effective couplings, Eqs. I|2.9|) . 1|2.10[1 . (|2.11|) and Eq. (|2.17|l . that for pertur- 
bative values for the tree-level coupling parameters the predicted results for SNR/ISB are very stable even for very 
large temperatures (in units of the regularization scale Ma) which is due to the slow logarithmic change with the 
temperature. As an illustration, consider for example the tree-level coupling parameters that satisfy the boundness 
condition Eq. A = 7 x 1(T 5 , A = 5 x 1(T 4 and A = -6 x 1CT 5 and = = 2. For these values 

of parameters Eq. I|2.6|l is satisfied and the one-loop equations for the effective masses predict ISB or SNR, along 
the <j) direction, for m 2 ^ > or m 2 < 0, respectively. Fig. 0] shows that the boundness condition also holds true 
for the effective temperature dependent couplings. Fig. shows the quantity A^(T), defined by Eq. (|2.16l) . which 
remains negative for the whole range of temperatures considered, thus predicting SNR/ISB along the <j) direction, in 
accordance with the WRG results 3]. At the same time, A^(T), shown in Fig. remains always positive. 

Note that the apparent almost constancy in a wide range of temperatures seem from the Figs. [S] and [S] is only 
a consequence of the effective couplings be only logarithmically dependent on T and the very small values for the 
tree-level couplings that we have considered. Had we taken larger values for the tree-level couplings, obviously would 
lead to a much larger variation with increasing temperature. 

Given the results shown above for the relativistic case, we can conclude, therefore, that the inclusion of thermal 
effects on the couplings does not exclude the possibility of SNR/ISB occurring at high temperatures. We recall that 
although the results were obtained with the one-loop approximation this feature does not seem to be an artifact of 
perturbation theory as confirmed by the results produced by nonperturbative methods, such as the Wilson Renor- 
malization Group procedure used in Ref. [3j , as well as the optimized perturbation theory used in Ref. |^| , where not 
only thermal corrections to the couplings are accounted for but also to the masses (like in the Schwinger-Dyson or 
gap equations for the masses). 
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FIG. 4: The boundness condition, Eq. \Y2.2t . expressed in terms of the effective couplings, is shown in the vertical axis as the 
ratio R(T) = A^(T)A,/,(T)/[3A(T)] 2 . The temperature is shown in units of the regularization scale Mo while N$ = N^p = 2. 
The tree-level values A^ = 7 X 10~ 5 , A^ = 5 x 1CP 4 and A = -6x 10 -5 were considered. 
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FIG. 5: The quantity A^(T) as a function of the temperature for the same values of parameters considered in Fig. 2] 

III. SEARCHING FOR SNR/ISB PATTERNS IN THE NONRELATIVISTIC CASE 

We now turn our attention to the analysis of similar SNR/ISB phenomena displayed by the relativistic model, given 
by the lagrangian density, Eq. I|2.1|l . in the case of its nonrelativistic counterpart. Let us first recall some fundamental 
differences between relativistic and nonrelativistic theories that will be important in our analysis. Firstly, the obvious 
reduction from Lorentz to Galilean invariance. Secondly, it should be noted that in the nonrelativistic description 
particle number is conserved and so, only complex fields are allowed. This second point will be particularly important 
to us since, for the processes entering in the effective couplings shown in Figs. n0 an dEl oni y those that do not change 
particle number (the elastic processes) will be allowed (e.g. this selects the processes (a) shown in Fig. |3]but not the 
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FIG. 6: The quantity A^(T) as a function of the temperature for the same values of parameters considered in Fig. 01 



(b) and (c), inelastic, ones). Another important difference between relativistic and nonrelativistic models concerns 
the structure of the respective propagators. While the relativistic propagator allows for both forward and backward 
particle propagation (which is associated to particles and anti-particles, respectively), the nonrelativistic propagator 
of scalar theories at T — only has forward propagation (see e.g. the discussion in Ref. 0]). Note however that the 
structure of the propagators (or two-point Green's function) in a thermal bath includes both backward and forward 
propagation |20|. which can be interpreted in terms of excitations to and from the thermal bath (or, equivalently, 
emission and absorption of particle to and from the thermal bath pljl. 

We should also say that, alternatively to the derivation of the nonrelativistic analog of i|2.1|) . we could as well 
consider the relevant equations leading e.g. to the derived effective couplings in the previous section and take the 
appropriate low-energy limit for those equations. However it is more practical, and indeed it is the procedure usually 
adopted in atomic and low energy nuclear physics, to start directly from the nonrelativistic Hamiltonian or Lagrangian 
densities. This is an one step procedure leading, say, to the Feynman rules that can be applied to any other quantity 
that we may be interested in computing, without having first to compute the corresponding relativistic expressions 
and then working out the corresponding nonrelativistic. So, let us now initially consider the nonrelativistic limit of 
the lagrangian density given by Eq. I|2.1(l . This can be obtained by first expressing the fields </> and tp in terms of 
(complex) nonrelativistic fields $ and VP as [l9l I22I |23| 



</>(£, t) 



1 



[exp(— imj,i)$(x, t) + exp^m^t)^* (x, t)] , 



and 



(3.1) 



ip(x, t) 



[cxp 



(— im^t)^(x, t) + exp(im^t)ty*(x, t)] 



(3.2) 



where it is assumed that the fields $ and \1/ oscillate in time much more slowly than exp(zm^t) and exp(zm^t), 
respectively. By substituting (|3.1() and l|3.2() in (|2.1() and taking the nonrelativistic limit of large masses, the oscillatory 
terms with frequencies and can be dropped. The resulting lagrangian density in terms of \1/ and complex 
conjugate fields becomes 



£($*,$, **,*) = _L \-i m Jd t ^*)^ + im^*(d t ^) - |V$| 2 + |ft$f 
2m^, L 



16m^ 



($*$) 
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1 



2m^ 
A 



-im^dt^*)^ + im^*(d t v) - |w| 2 + \d t ^\ 2 ] - y^rC***) 2 

, (3.3) 



where we have assumed for simplicity, in the derivation of the last term in (|3.3p . the cross- fields interaction term, that 
rrirf, ^ to,/, 3 . By further considering 

|d f $| 2 « 2m Im($9 t $*) , 

|9 t *| 2 < 2m^Im(*a t $*) , (3.4) 

we can omit the terms with two time derivatives in Eq. (|3.3|l . So the Lorentz invariance in (|3.3|) is lost and the 
nonrelativistic analogue of (|2.1[1 is obtained. The interaction terms in Eq. (|3.3|l are the same as those obtained by 
approximating the usual nonrelativistic two-body interaction potentials by hard core (delta) potentials, e.g., 

J d 3 x$*(x,t)$(x,t)V^,(x-x')$*(x',<)$(x',t) [$*(x,t)$(x,i)] 2 , 

J d 3 a;**(x,i)*(x,i)V*(x-x')**(x',i)*(x',i) 5 * [^*(x, <)*(x, t)] 2 , 

J d 3 x$*(x,t)$(x,i)t/ <E >v[,(x-x')4'*(x',<)*(x',i) -> g^ [$*(x, t)$(x, t)] [tf*(x,t)tf (x, t)] . (3.5) 

The approximation of the two-body potential interactions like in (|3.5J) is also commonly adopted in the description 
of cold dilute atomic systems, where only binary type interactions at low energy are relevant. In that case, the local 
coupling parameters g$, and </<j>$ are also associated to the s-wave scattering lengths [jj, e.g., gi = 2-Kai/mi. For 
nonrelativistic systems in general, besides the two-body interaction terms like 1)3.5(1 (in the hard core approximation) 
we can also include additional one-body like interaction terms, e.g., etc. This is the case when we submit the 

system to an external potential (for example a magnetic field). It can also represent an internal energy term (like the 
internal molecular energy relative to free atoms in which case the fields in the lagrangian would be related to molecular 
dimers) . In models of superconductivity a constant one-body like interaction term represents the opening of an explicit 
gap of energy in the system. In the grand-canonical formulation can represent chemical potentials included in the 
action formulation, so that one can also describe density effects (in addition to those from the temperature). In order 
to retain the symmetry breaking analogies to the previous relativistic model l|2.1[> . and since our intention here is to 
keep the analysis as general as possible, we shall also consider additional one-body interaction terms for <I> and \t 
while their precise interpretation is left as open and will depend on the particular system Eq. (|3.3|) is intented to 
represent. 

With the considerations assumed above, we therefore take the following nonrelativistic lagrangian model, that is 
analogue to the relativistic model Eq. (|2.ip . 



£($*, $,**,*) = $* [ id t + — !— V 2 ^($*$) 2 

2m$ / 3! 

1 T-,2 \ , T , ,,_,T,*,T, 5* /lTr *, T ,N2 



+ id t + -z V 2 * - k^*^ - Trf(***) 

\ J 3! 

- £/($*$)(***) , (3.6) 

where we have expressed the derivative terms in their more common form (by doing an integration by parts in the 
action context). The numerical factors and signs in the one and two-body potential terms in l|3.6|) have been chosen in 
such a way so that the potential in l|3.6|) are analogous the one considered in l|2.1ll . The coupling constants shown in 
(|3.6() are related to those in l|3.3(l by gi = 3Xi/(8mf) (with i = Vt) and g = A/(4rn$m*) while m$ and to* represent 



3 Note that for equal masses there is the possibility of an additional symmetric interaction term of the form [(<E>'J'*) 2 + (<t'*>I') 2 ] in 13.31 . 
however this term will not be relevant for our analysis and conclusions since it can be absorbed in a redefinition of the cross-coupling 
constant A, especially when we work with densities, or averages of the fields, like in an effective potential calculation. 
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the (atomic) masses. In addition, notice that for the nonrelativistic limit which leads to Eq. I|3.6[l to be valid, one 
must keep T <C rrii. Since for nonrelativistic systems in general, the masses rrii are of order of typical atomic masses, 
rrii ~ 0(1 — 100)GeV, and the typical temperatures in condensed matter systems are at most of order of a few eV, 
this condition will always hold for the ranges of temperature we will be interested in below. 

For multi-component fields, Eq. I|3.6[) is the nonrelativistic multi-scalar model with symmetry U(N<s>) x U(N^) that 
is the analogue of the original relativistic model Eq. p.ljl . For simplicity, in the following we assume the simplest 
version of (|3.6[1 where iV$ = = 1, corresponding to an U(l) x U{1) symmetric model. In this case, by writing the 
complex fields in terms of real components, 



* = A= (0i + ifo) , * = 4= W>i + #2) , (3.7) 

we see that the lagrangian model (|3.6[) falls in the same class of universality as that of Eq. (|2.1[1 for the case of 
the 0(2) x 0(2) symmetry. The extension to higher symmetries can be done starting from l|3.6[) but the case of 
simplest symmetry involving the coupling of complex scalar fields will already be sufficient for our study (physically, 
this system may, for example, describe the coupling of Bose atoms or molecules in an atomic dilute gas system) . 

Just like in the relativistic case, in the current application we consider the Ki, that appears in Eq. (|3.6|) . as simple 
temperature independent parameters for which thermal corrections arise from the evaluation of the corresponding 
field self-energies. Now, to make contact with the analogous potential used in the prototype relativistic models for 
SNR/ISB, we take the overall potential as being repulsive, bounded from below. This requirement imposes a constraint 
condition analogous to the one found in the relativistic case, g-q, > 0, gq, > and g^g® > 9g 2 . At a given temperature, 
the phase structure of the model is then given by the sign of Ki{T) = /c, + Sj, where £j is the field temperature 
dependent self-energy. In the broken phase Ki(T) < 0, while in the symmetric phase Ki(T) > 0. The phase transition 
occurs at Ki(T = T£) = 0. At the one-loop level the diagrams contributing to the field self-energies are the same as 
those in the relativistic case. The main difference is that the momentum integrals in the loops are now given in terms 
of the nonrelativistic propagators for <E> and 'J, 



AK,q) = — , > , (3.8) 

-iu„ + Wi(q) 

where w„ = 2%nT (n — 0,±1,±2, . . .) are the bosonic Matsubara frequencies and Wi(q) = q 2 /(2mj) + «j. At the 
one-loop level we then obtain (see appendix) 



« — — ran x - 



n— — 00 



2ff*/3 | 9 
-iu„ + Wi(q) -iuj„ + Wj(q) 



(3.9) 



The sum in Eq. 13. 9|) can be easily performed and the resulting momentum integrals lead to well known Bose integrals 
(see e.g. Ref. |2^| and also the appendix for the derivations of these equations). We then obtain the results 



K$(T) = k$ + [ ) Li 3/2 [exp(-K$/r)] +g ( ) Li 3/2 [cxp(-K*/T)] , 



3 V 27r / 



V 2tt J 



(3.10) 



and 



K * (T) = K<i H ^— f ^ J Li 3/2 [exp(-K*/T)] +g f j Li 3/2 [exp(-K$/T)] , (3.11) 

where Li„(zj) is the polylogarithmic function. Like in the relativistic effective mass terms, in Eqs. (I3.10|) and (I3.11|) 
we have once again limited to showing the temperature dependent corrections coming from S^, omitting the divergent 
(zero-point energy terms) contributions to the effective one-body terms (so the tree-level parameters in Eqs. (|3.1U|) and 
(|3.11l) are assumed to be already the renormalized ones). Let us now consider the high temperature approximation, 
which in the nonrelativistic case is valid as long as re, <C T <C m,. Then, one is allowed to consider the approximation 
for the polylogarithmic functions in Eqs. (|3.1U|) and 13.11|) . Li 3 / 2 [exp(— Ki/T)] ~ C(3/2), where ((x) is the Riemann 
zeta function and C(3/2) — 2.6124. One can then write the two equations (|3.1U|I and l|3.11|l in the high temperature 
approximation more compactly, as 
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T 



3/2 



Ki (T) ~ + C(3/2)Af R , (3.12) 

where we defined the quantity Af R analogous to that of the relativistic case, 

^=\ 9i mr+9mT , (3-13) 

in terms of which we obtain the critical temperature for symmetry restoration/breaking analogous to the relativistic 
expression, 



Af R C(3/2) 



2/3 

(3.14) 



Eq. (|3.14|l shows that there are three interesting cases which depend on the sign and magnitude of the cross coupling 
g. Taking < and g > one observes a shift in the critical temperatures indicating that the transition occurs at 
lower temperatures compared to the decoupled case (g = 0). If g < but \g\ < (2/3)gi(mi/rrij) 3 ^ 2 then the transition 
occurs at higher temperatures. Despite these quantitative differences symmetry restoration does take place in both 
cases. Now consider for instance the case, with k$ < 0, where g < but \g\ > (2/3)(7$(m$/m^) 3 / 2 (in this case, and 
assuming m$ ~ m*, the boundness condition assures that \g\ < (2/3)g*). Under these conditions, for the $ field 
we have a similar situation as that for the corresponding relativistic case studied in Sec. II, where, Eq. (|3.14|l does 
not give a finite, positive real quantity. This is a manifestation of ISB (for k$ < 0, or SNR, for k$ > 0) within our 
two-field complex nonrelativistic model being analogous to what is seen in the relativistic case. At the same time the 
field 4" suffers the expected phase transition at a higher T c compared to the g — case. As in the relativistic case, 
which field will suffer SNR/ISB depends on our initial choice of parameters in the tree level potential, so it is model 
dependent. The SNR/ISB result as seen above in this nonrelativistic model is however misleading, as we next show 
by considering the same phenomenon in terms of the effective, temperature dependent, coupling parameters. 

Physically, the possibility of SNR/ISB occurring at high temperatures as predicted by the naive perturbative 
approximation, Eq. (|3.12|) . is much harder to be accepted in the nonrelativistic case (condensed matter) than in the 
relativistic case (cosmology) , which already seems to indicate that the results here should be different when the simple 
perturbative calculations are improved. 

Following the analogy with the relativistic model calculation done in Sec. II, we next evaluate the temperature 
effects on the nonrelativistic couplings. 

The diagrams contributing to the effective self-couplings g$(T) and g^,{T), at the one-loop level, are those shown in 
Figs, n an d 121 respectively, except that the s-channel ones with internal propagators for fields different from those in 
the external legs (the second one-loop diagrams in Figs. ^ and [21 made of vertices nonconserving particles) are absent. 
For the effective cross-coupling g(T), as discussed at the beginning of this section, the diagrams contributing at the 
one-loop level are the ones shown in Fig. IHIa), corresponding to the particle number preserving processes. Using i|3.8|) 
for the nonrelativistic field propagators, the explicit expressions for the effective couplings at one-loop order are found 
to be (see also the appendix) 



g 2 f d 3 q 1 

9t>{T) = 9<s> ~ -y J (2?r)3 2uj>(q) ^ + 2n(yUJ ^ + n ( UJ <s>) i 1 + n {v<s>)]} 

- 3g 2 J ^L/3n(o;*)[l + n(^)] , (3.15) 

- 3g 2 J A^/3n(w $ )[l + n(^)] , (3.16) 



and 
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2.9 f d 3 q 

g( T )=9 - y J T^Xs/ 3 {fl 1 *"^*) I 1 + n(uj< s ,)} + gyn(ujit) [1 + n(w t )]} 

The numerical factors in Eqs. 1)3.15(1 . 1)3.16(1 and 1)3.17(1 are due to the symmetries of the diagrams and normalizations 
chosen in Eq. 1)3.6(1 . n(uj) in the above equations is the Bose-Einstcin distribution. 

The zero temperature contributions in Eqs. ((3.15)1 - l|3.17|l are divergent and require proper renormalization. This 
is mostly simply done by performing the momentum integrals in d = 3 — e dimensions and the resulting integrals are 
all found to be finite in dimensional regularization (when taking e — > at the end) . The momentum integrals for the 
finite temperature contributions in Eqs. 1)3.15)1 - 1)3.17)1 can again be performed with the help of the Bose integrals 
given in the appendix. We can further simplify the equations by considering parameters such as to$ ~ = to, 
and temperatures satisfying k < T < m in which case the zero temperature corrections to the couplings 
are negligible compared to the finite temperature ones and can safely be neglected. At leading order, in T, we then 
obtain the results 

g*(T) ~ <?$ - — ^ — (5 9 | + 9g 2 ) + 0{k/T) , (3.18) 
g*(T) ~ .9* - (5gl + 9g 2 ) + 0(n/T) , (3.19) 

12-7T V K 



l( T ) —9 - -tt \ —-9 (g + -ir + ) + 0{k,T) . (3.20) 



and 



4tt V k V 3 3 

Note from Eqs. (|3.18() . (|3.19(l and l|3.20|l that the effective couplings in the nonrelativistic theory have a much stronger 
dependence with the temperature than in those in the equivalent relativistic theory, Eqs. 1)2.9)1 . 1)2.10)1 and 1)2.11)1 . We 
therefore expect to see larger deviations at high temperatures for the effective couplings as compared with the same case 
in the relativistic problem (by high temperature we mean here temperatures larger than the typical one-body potential 
coefficients in 1)3. 6(1 . but much less than the particle masses, see above). It is also evident from the analysis of higher loop 
corrections to the effective couplings in the nonrelativistic model that all bubble like corrections contribute with the 
same power in temperature as the one-loop terms, which can easily be checked by simple power-counting in the momen- 
tum. A side effect of this is the breakdown, at high temperatures, of the simple one-loop perturbation theory applied 
here. Another symptom is the apparent running of the effective self-couplings, shown above, to negative values for suffi- 
ciently high temperatures given by T > Tb roa k ~ min (l2ir^J n/ (2m 3 )g<s, / (5<?| + 9.g 2 ), YItv^/kJ (2m?) ,g*/(5.g| + 9g 2 )j . 

Nevertheless, it is easy to check that (for the parameters adopted below) the results obtained by just plugging Eqs. 
1)3.18)1 . 1)3.19(1 and 1)3.20)1 above into Eqs. 1)3.12(1 already show a drastic qualitative difference between this simple im- 
proved approximation and the naive perturbative evaluation given by Eq. ((3.12)) . It looks that for the nonrelativistic 
case SNR/ISB is a mere artifact of perturbation theory. Intuitively, this is already a rather satisfactory result which, 
as we will show below, will be confirmed by a nonperturbative resummation of the bubble like corrections. The re- 
summing of all leading order bubble corrections to the couplings can again be done by solving the set of homogeneous 
linear equations for g$>(T), gv(T) and g(T), 

9<s>( T ) = 9<s> - 9<s>( T ) y h(@K<!>) ~ 3g(T)g > 
= 5*-5*( T )y hip*,*) - 3g(T) g I 2 (/3n$) , 

g(T) = g - g{T) ^ I 2 (/3k $ ) - g{T) ^ I 2 (/3k*) - g 9 (T) ^ h(P^) - ff*(T) ^ hip**) , 

- g(T)gI 3 (p K ^p K9 ) , (3.21) 

where we have defined the functions 
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(2tt)3 2wi(q) 
d 3 q 



0n(u>i) [1 + n(tJi)] 




d 3 q 1 

{1 + 2n(uji) + %Pu)iTi{wi) [1 + n(uji)}} , 



1 + 2n(wi)] - 2 _ 2 [1 + 2n(o; i )] ^ . (3.22) 



One is now in position to investigate how thermal effects on the effective nonrelativistic couplings manifest 
themselves in phenomena similar to SNR/ISB. First we show in Fig. [7| some representative results for the ef- 
fective interactions obtained from the solutions of Eq. (|3.21|l . The tree- level parameters considered here are: 
S$ = 2x lCT^eV- 2 ,^ = lCT 16 eV~ 2 ,.g = -10- 16 eV~ 2 , m$ ~ = lGeV and k$ = rc w = IneV. We note 
that all couplings tend to evolve to zero at very high temperatures as T gets closer to m. So, the apparent instability 
caused by the fact that g^(T) as well as g$>(T) could become negative beyond some temperature, Tb re ak, as suggested 
by Eqs. (|3.18() and (|3.19|l . has disappeared completely when the nonperturbative flow of the couplings are considered 
(using the tree-level parameters given above and from Eqs. H3.18|l and (|3.19(1 . one would get Tb rca k ~ 0.019 eV (or 
~ 220 K)). 

It is tempting, by looking at Fig. to associate this to a free model at high energies, however we must recall that 
the nonrelativistic model, Eq. (|3.6[1 . will eventually no longer be valid for such high values of temperature, in which 
case it should be replaced by the original relativistic model (in any case the model Eq. i|3.6[l should of course be 
regarded as an effective model valid at low energy scales only). 
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FIG. 7: The non relativistic effective couplings as a function of temperature (shown in units of eV). 

In Fig. we show the equivalent of Eq. i|XT3|) for the $ and * fields, A| R (T) and A^ R (T), respectively, given in 
terms of the temperature dependent bubble resummed couplings for the same tree-level parameters considered above. 
We note that for the parameters considered A^ R (T) is initially negative and reverse sign at some temperature, 
indicating that symmetry breaking at high temperatures tend to happen in the W field direction, while the potential 
in the $ direction remains unbroken (actually, the temperature where A^ R (T) crosses zero, is close to the point of 
a reentrant transition for 'J). A^ R (T) however always remain positive, which then points to no transition in the $ 
direction. These aspects are also clearly seen in the plot, shown in Fig. ^\ for the effective one-body term k^(T), 
expressed in terms of the effective nonperturbative couplings. In terms of the parameters considered, symmetry 
breaking is seen to happen at a temperature T^ S ^ ~ 3.4 x 10" 4 eV (or ~4K) while the reentrant phase (symmetry 
restoration) happens at a temperature T^jg ~ 1.4 x 10~ 2 eV (or ~ 161 K). In between these two temperatures we 
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see a manifestation of an ISB phase. In the $ direction there is no symmetry breaking or reentrant phases at any 
temperature for the parameters considered. 
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FIG. 8: The quantities A^ R (T)/m 3 ' 2 and A^ R (T)/m 3 ' /2 as a function of the temperature for the same values of parameters 
used in Fig. [7] 




T(eV) 

FIG. 9: The nonperturbative effective one-body term k*(T) as a function of the temperature for the same values of parameters 
used in Fig. Both quantities are expressed in units of eV. The arrows indicate the points of SB and SR, with an intermediary 
ISB phase happening between the temperatures T^ B < T < T C SR . 

In Fig. 1101 we show a phase diagram for the system as a function of the tree-level coupling g^ and the temperature. 
The thin horizontal line at = 10~ 16 eV -2 illustrates the reentrant transition (through an inverse symmetry breaking) 
shown in Fig. All other parameters are the same as considered above. Note that the condition of stability, 
g^g$> > 9g 2 , expressed in terms of the nonperturbative and temperature dependent couplings, is always satisfied at 



Hi 



any temperature for the parameters considered for the previous figures (the effective couplings and 5$ also remain 
always positive, as is clear from Fig. [7|and previous discussion). 
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FIG. 10: A phase diagram of the system for fixed parameters: g$ — 2 x 10 15 eV 2 ,g = —10 16 eV 2 , mg> ~ = 1 GeV 
and k<i. = Kg< = IneV. The dotted line separates the region of stability (R NR (T) > 1) from that of instability (R NR (T) < 1), 
according to the value of the ratio R NR (T) = gy(T)g<s>(T)/[9g 2 (T)]. The regions of SR and SB, in the \& direction, are also 
shown. The thin horizontal line indicates the value of gm used in the text. 



IV. CONCLUSIONS 



We have reviewed how symmetry non-restoration and inverse symmetry breaking may take place, at arbitrarily large 
temperatures, in multi-field scalar relativistic and nonrelativistic theories. These counter intuitive phenomena appear 
due to the fact that the crossed interaction can be negative while the models are still bounded from below. We have 
recalled that, in the relativistic case, SNR/ISB are not a mere artifact of calculational approximations. We have then 
set to investigate the possible SNR/ISB manifestation and consequences in a nonrelativistic U(N$) x U(N^) scalar 
model of hard core spheres by considering the simplest case, JV$ = N-q, — 1, which may be relevant for condensed 
matter systems of bosonic atoms or molecules. 

Performing a naive perturbative one-loop calculation, which includes only the first thermal contribution to the self- 
energy, we have shown that, for negative values of the crossed coupling, SNR/ISB can take place like in the relativistic 
case. However, the manifestation of SNR/ISB in condensed matter systems of hard core spheres seems to be more 
counter intuitive than in the relativistic case where the model may represent, e.g., the Higgs sector. With this in 
mind we have investigated the explicit (temperature) running of the nonrelativistic couplings. One first improvement 
was to evaluate the perturbative one- loop thermal corrections to the couplings which already indicate that SNR/ISB 
do not seem to happen, at high temperatures, for the nonrelativistic case. Next, we have resummed the bubble like 
contributions in a nonperturbative fashion. This procedure fixed the instability problem related to the possibility of 
gi(T) becoming negative as observed in the calculation which considered only the simplest one-loop corrections to the 
couplings. Our nonperturbative calculation also showed that the phase transitions happening in the nonrelativistic 
case includes a continuous SB/SR pattern characterized, at intermediate temperatures, by a reentrant, continuous 
transition. Therefore, we can state as our major result that, contrary to the relativistic case, SNR/ISB does not seem 
to occur in the nonrelativistic model of hard core spheres. Instead, reentrant like phenomena become possible, as our 
results have indicated. 

Finally, it would be interesting to investigate SNR/ISB in connection with the Bose-Einstein condensation problem 
and the present work gives some of the ideas and tools needed for this task that we hope to pursue and report in a 
future publication. 
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APPENDIX A: TEMPERATURE DEPENDENT CORRECTIONS FOR THE NONRELATIVISTIC 

MODEL 

Here we give the main steps used in the evaluation of the thermal masses and couplings for the nonrelativistic 
model. Similar derivations for the relativistic model can be found e.g. in Refs. 1 11 1 1 . 

The effective one-body terms ftj(T) and the effective couplings (two-body terms) for the nonrelativistic model are 
most easily to obtain directly from a computation of the one-loop effective potential for the model lagrangian (|3.tjfl . 
The effective one and two-body terms will then be identified with the appropriate derivatives of the effective potential. 

As usual in the computation of the one-loop potential, we start by decomposing the fields $ and '5 in H3.6|l in terms 
of (constant) background fields (which, without loss of generality can be taken as real fields) </>o and i/io, respectively, 
and fluctuations <f> and ip, which in terms of real components, become 



v2 



(Al) 
(A2) 



When substituting Eqs. (|A1|) and (| A2|) in (|3.t)[l we only need keep the quadratic terms in the fluctuation fields for the 
computation of the one-loop potential for the background fields 0o and ipa- We then obtain the (Euclidean) lagrangian 
density in terms of <f>o and ipo, 



K$ 2 yq, 4 n, v 2 yw 4 y 2 



.9* 



9 



4! 



4! 



1 



-\ ■ M ■ x + cubic and quartic interaction terms 



(A3) 



where we have defined the vector \ 
fluctuations, 



h\, (f>2, tpi, ^2) and M is the matrix operator for the quadratic terms in the 



I 2m$ 



M = 



V 



-id T 

Mo 




2m$ 



id. 

K <S> T — Vo 






9* ^ ^ £^,2 



-V 

2 m a- 



(A4) 

The partial time derivative in l)A4|l is over Euclidean time: d T — d/dr, t — it. By performing the functional 
integration in the quadratic fluctuations \i the one-loop effective potential V e ff(4'Oi' t l'o) obtained from Eq. (|A3|) is 
given by 



g<j>oi>o 



— i$ r 
(9/<9t, r 



2vn q> 








9* ,/,2 4. 9^2 
g ^0+2^0 



VcS (fc, = f ^ + § ^ + + f ^ + + 5 ^ det M , 

where the last term on the rhs of (|A5fl comes from the functional integral over the components of x, 



(A5) 



1 



In det M 



1 



DcfriDfaDipiDfa exp 







<r,- ( I x . m . x 



(A6) 



and V is the volume of space. Expressing Eqs. (|A6|) and (|A4|) in the space-time momentum Fourier transform form, 
we then also obtain 
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In det M 



9 



d 3 q 

(gipl [u; # (q) - 



2 



ln{[w*+i£(q)] [c4 + £|(q) 



9<s> 



w*(q) 



5* ,2 i 9 ,2 



with 



#*(q) = 



and 



23*(q) = ^ |^(q) + 9 f^ 2 + |$j [^(q) + ^ 2 + |$ 
From V^ff(</>o,^o)) Eq. (|A5|) . we now define the effective one-body terms as 
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K$(T) 



««r(T) = 



d 2 V c $(<j)o,il>o) 



d 2 V cB (cl)Q,ip ) 



dip 2 



00=0,^0=0 



which then gives 
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~3~P^J (27r)3^+c4(q) ' 



E 



<i 3 q w$(q) 



4?/ 



and 



/^s 2 5>f 1 ^ /■ d J ( 



d 3 q w*(q) 1 



d 3 q w*(q) 
(27r)3^ + w |(q) 



d 3 q w$(q) 



(A7) 

(A8) 

(A9) 

(A10) 
(All) 

(A12) 
(A13) 



The Eqs. jUg) and (jA"l3|) can also easily be expressed in terms of the free nonrelativistic propagators D$(w„, q) and 
Dy,(uj n ,q). The sum over the Matsubara frequencies in (|A12|) and i|A13(l are easily performed by using the identity 



1 ^ 1 1 

B Jl + uj 2 u) 

1 il 



(A14) 



where n(w) is the Bose-Einstein distribution, 



n(u>) 



(A15) 



e?" - 1 ' 

The momentum integrals in Eqs. I|A12|I and (|A13I) can be expressed in terms of standard Bose integrals as follows 
(see for example j24|). Consider the integral (where uj — q 2 / (2m) + k and 77 = (3k) 
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(A16) 
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where we used the definition for the polylogarithmic function, 



(A17) 



;=i 



Another useful momentum integral that also can be obtained from <|A16|) is 



/ 



d 3 q 
(27T) 



/ m \ 3/2 
T n{u)[l + n(u>)] = ( ^3 J Li i/2 (e 



(A18) 



We also have the results obtained from the the polylogarithmic functions in the high temperature approximation, 
k <C T, and that are used in the text, 



Li 3/2 {e- 0K ) = C(3/2) - 



(A19) 



and 



T 



Lii/2 (e-"") = \U- - C(l/2) + C(«/T) . 



(A20) 



Using (|A16|I in Eqs. (|A12|) and (|A13|) we obtain the results quoted in the text, Eqs. I|3.10j) and l|H.llj) . 

The two-body effective terms are also defined analogously as 
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and 



9 4 y off (0 o , 1P0) 
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(A23) 



From Eqs. I|A5[) and l|A7|l we then obtain for Eqs. (|A21I) - l|A23|l the results 



d\ J 2 1 8 2 c4(q) 

(2tt)3 1 5 *u£ + <4(q) 3 ff %2 +c4(q)] 



3^ 



<4(q) 



6 5 2 



<4(q) 



•<4(q)r 



(A24) 



= .9* + 



d 3 c 



3 5 2 - 



(2tt)3 l 3 *^+^|(q) 3 



8 



5*- 



<4(q) 



■<4(q) 



65 



<4(q) 



+ <4(q) [<4 + c4(q)r 



(A25) 



and 



20 



1 ^ f d 3 q J 2 1 4_ w|(q) 



9{T) = 9 + ~ V / 77^ < o993> ,. 2 , -tt-x " o55*" 



/3^7 ( 27 r)3 |3^^ +w | (q) 3 ^ K+W | (q)]2 

2 1 4 cj|(q) 2 cj $ (q)cj^(q) 

3 5ff *^ + w a (q) 3 3S* K + w|(q)]2 3 [W 2 + w a (q)] ^ + w a (q)] 



The Eqs. (|A24fl - (|A26(1 again can be expressed in terms of the free nonrelativistic propagators D$(u n ,q) and 
D<p(o; n ,q), which can then be identified with the corresponding one-loop diagrams that contribute here, depicted in 
Figs 1, 2 and 3a. All sums over the Matsubara frequencies in (|A24|) - (|A26|I are again evaluated with the help of the 
identity 1A14I) . from which we obtain the results shown in Eqs. (|3.15|) . I|3.16[) and l|3.17|) . 
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